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1.1.  Introduction 


Let  us  consider  a  continuous,  irrotational,  two-dimensional  flow 
limited  by  a  free  surface.  The  fluid  is  supposed  to  be  incompressible, 
perfect  and  heavy.  The  complex  potential  of  the  flow  can  be  written 
as  follows  : 

F(z)  »  Voz  +  f(z)  “ 

where  '^o  represents  the  velocity  of  the  uniform  flow  at  upstream 
infinity  amd  f(z)  the  complev  perturbation  potential* 

We  intend  to  consider  the  stream  function  : 


y  .  Voy  .  y 

Where  ^  represents  the  stream  perturbation  which  is,  as  ,  an 
harmonic  function.  This  means  that  the  ^  potential  meets  the 


equation  : 

2  2 
^  ^  ^  *  0 

^  ^  y^ 


12,  Generalities  on  the linear  theory  of  the  small  perturbations 


for  a  flow  limited  by  its  free  surface. 


Let  us  consider  a  two-dimensional  flow,  as  defined  above,  where  a 
thin  profile  represented  by  its  chord  AF  =  S  is  placed  at  a  distance  £ 
of  the  static  level  of  the  free  surface  (fig.  1.1.). 
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The  free  surface  is  a  stream  line  of  the  flow  (normal  velocity 
equal  to  zero)i  The  Bernoulli  theorem  applied  to  this  linei  between 
the  two  sections  -  oo  and  +  ,  gives  ; 

po  +  Vo^  “  £  +  gy 

P  2  P  2 

where  po,  Vo  are  the  pressure  and  the  flow  velocity  at  upstream 
infinity,  and  p  and  V  the  same  elements  in  the  x  section  of  the 
perturbed  flow.  In  the  expression  (l ,4.)  we  insert  the  components 
u,  V  of  the  perturbation  velocity. 

^  +  J/o  =  £  +  Vo^  2Vou  u^  »  v^  +gy  (l»5») 

P  2  P  2 

Whence,  neglecting  the  second  order  terms  ; 

p  _  po  =  -  ^Vou  -  p  gy  i 

at  the  free  surface,  the  pressure  is  equal  to  po,  whence  the  condition  ; 

Vou=-gy  (l'7.) 

The  sliding  condition  on  the  profile  and  at  the  free  surface  shows 
that,  there,  the  normal  velocities  are  equal  to  zero,  this  being 
in  accordance  with  the  physical  reality  (always  with  respect  to  the 
linear  theory).  It  is  written  : 


dr: 


V 

Vo  +  u 


V 

Vo 


(1.8.) 
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We  shall  insert  to  1.7.  and  1,8.  the  perturbation  stream  function, 
by  meauas  of  already  known  relations  : 


So  : 


and, 


'd  y 


V  a  -  ^  ^  (1.S.) 

d  X 


Vo  ^  ^  »  -  gy 

d  y 

dV  a  .  1  , 

dx  ^  X  Vo 


(1 .10.) 


(1.11.) 


Taking  this  last  condition  projected  on  to  the  x's  axis  and  integra¬ 
ting  between  -  oo  and  x,  we  obtain  : 


y 


si 


Vo 


(1.12.) 


where  will  represent  the  ordinate  at  the  free  surface.  If  we 
insert  (l,12.)  in  (l.lO.)  ; 


a  g  ,  (1.13.) 

d  y  vo^ 

The  sliding  condition  on  the  linearized  profile  in  the  x's  direction, 
integrated  between  the  leading  edge  and  the  abscissa  ,  allows  us 
to  obtain  : 

y  ( I  )  -  -  A  (1.14.) 

Vo 
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To  the  assumptions  which  have  enabled  us  to  define  the  (l.12.), 

(l.13.),  (l«14.)  conditions,  we  add  the  assumption  of  the  linearization 
of  the  boundary  conditions  ;  these  conditions  are  written  in  projec¬ 
tion  and  fulfilled  on  the  x's  axis.  The  condition  (l.13.)  will  now  be 
written  ; 

(1.15.) 

,  — V 

(  n  normal  directed  towards  the 
inside  of  the  fluid) 


0 't' 


d  n 


Vo 


and,  by  inserting  the  factor 


we  shall  write  ; 

-  =  r  Y 


(1.16.) 


(1.17.) 


We  still  have  to  add  another  condition  which  takes  into  account  the 
characteristics  of  the  flow  at  upstream  of  the  profile.  In  this  region, 
there  are  no  waves,  while  they  may  appear  at  downstream  of  the  profile. 


1.3.  Set-up  of  the  problem. 

The  solution  of  the  direct  problem,  i.e.  the  study  of  a  profile 
whose  characteristics  in  an  Indefinite  flow  are  known,  when  it  is 
located  in  a  flow  limited  by  a  free  surface,  has  been  envisaged  by 
various  authors, 

Vladimirov  (l)  quotes  various  works  permitting  to  determine  the  total 
efforts  acting  on  an  immerged  profile.  The  author  himself  and  the 
people  mentioned  in  his  article  do  not  take  interest  in  the  details 
of  the  flow  -  distribution  of  pressures  and  velocities  -  though  this 
problem  be  of  a  great  importauice  for  the  study  of  all  the  cavitation 
phonoMna . 
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More  recent  experimental  and  theoretical  researches  have  tried  to 
clear  up  this  problem.  Among  the  former,  let  us  mention  the  works  of 
Ausman  (3)  and  Parkin,  Perry  and  Wu  (4).  From  a  theoretical  point 
of  view,  the  work  of  Nishiyamn  (2)  and  of  Isay  (5)  presents  a  great 
interest. 

It  is  obvious  that  any  solution  allowing  the  determination  of  the 
shape  of  the  profile  for  a  distribution  of  imposed  pressure  (inverse 
problem)  should  be  welcome  by  the  engineers  in  charge  of  its  design. 

Two  works,  Nishiyama  (6)  and  Parkin  and  Peebles  (7)  deal  with  this 
problem.  But  the  former  only  taJces  into  account  the  free  surface  effect* 
Parkin  and  Peebles  start  from  the  assumption  that  the  profiles  in 
operation  are  placed  at  a  depth  which  is  equal  to  twice  their  chord 
value,  there  where  the  surface  effects  can  be  neglected  without 
Introducing  important  errors  into  the  computation. 

The  utilization  of  mathematic  models,  well  known  in  aerodynamics  in 
order  to  solve  these  problems,  has  logically  led  us  to  consider  the 
possibility  of  using  the  analog  computation  methods  -  which  have 
already  proved  satisfactory  in  aerodynamics  -  to  find  out  the  nume¬ 
rical  solution  of  this  problem.  We  thus  avoid  a  great  number  of 
unsuccessful  attempts,  especially  while  dealing  with  the  inverse 
problem.  We  consider  that  the  flow  generated  by  a  profile  is  due 
to  the  superposition  of  two  effects,  the  lift  effect  and  the  thickness 
effect.  The  former  is  equivalent  to  a  vortex  distribution  along  the 
projection  of  the  profile,  according  to  the  x's  axis.  The  later  corres¬ 
ponds  to  a  distribution  of  sources  and  sinks  on  the  same  projection. 

The  intensity  of  the  vortices  is  directly  in  relation  with  the  local 
load, while  the  mass  flow  of  the  sources  and  sinks  is  related  to  the 
thickness  law.  The  local  load  and  the  thickness  law  arc  the  basis  of 
the  problem. 
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From  the  analog  point  of  view,  it  appeared  that  the  condition  at  the 
free  surface  was  very  difficult  to  represent .  We  thus  have  been  obli¬ 
ged  to  modify  the  approach  of  the  problem  in  order  to  obtain  an  easier 
representation,  by  converting  the  Fourier  condition  (1,17*)  which, 
under  this  form,  implies  the  use  of  negative  resistances,  into  a 
Dirichlet  condition  satisfied  by  the  superposition  of  solutions  ob¬ 
tained  auaalytieally  for  a  vortex  or  a  source.  These  free  surfaces 
are  such  that  they  always  meet  the  conditions  (l,12.)  and  (1.17.)» 

The  computation  process  that  we  are  going  to  adopt  in  order  to  solve 
the  inverse  problem,  presents  itself  as  follows  : 

a)  Computation  of  the  free  surface  for  the  chosen  vortices  distribu¬ 
tion  and  for  a  given  flow. 

b)  Determination  of  the  mean  line  of  the  profile  before  the  vortex 
distribution  and  the  free  surface  of  (a)  par. 

c)  Analog  determination  of  the  free  surface  for  the  chosen  sink*  and 
sources  distribution  and  for  a  given  flow. 

d)  Analog  determination  of  the  shape  of  the  profile  which  corresponds 
to  (c)  par. 

This  settles  down  the  plan  to  follow  in  the  course  of  the  next  paragraphs. 


2.  Free  Surface  of  .the  Perturbed  Flows, 


2.1.  Set-up  of  the  problem. 

In  thislsection,  we  shall  determine  the  shape  of  the  free  surface 
generated  by  a  perturbation  -  vortex  or  source  -  immerged  in  a  uniform 
flow.  The  calculation  is  effected  by  the  "image"  method  (8)  where"one 
suppose  the  existence  of  a  S  system  of  sources  and  sinks  in  a  fluid 
which  possesses  one  or  several  boundaries  C,  If,  when  placing  a 
S'  system  of  sources  or  sinks  in  the  external  region  of  C  and  permit¬ 
ting  the  incoming  of  the  fluid  in  the  whole  region,  we  find  out  that 
C  is  a  stream  line  ;  the  S'  system  is  called  the  image  in  C  of  the 
S  system". 
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In  our  particular  case,  the  free  surface  is  a  stream  line  limiting 
a  flow  where  is  placed  a  vortex  or  a  source  at  a  distance  -  f  from 
this  surface  and  at  x  =  0.  The  vortex  or  image  source  will  then  be 
located  at  a  distance  +  f  from  this  surface  and  at  x  =  0,  To  the 
equations  system  thus  obtained  we  shall  have  to  add  a  certain  function 
permitting  to  fulfill  the  condition  (1»17.)  at  the  free  surface. 

2.2.  Free  Surface  for  an  Immerged  Vortex. 

The  complex  potential  of  a  vortex  is  written  : 

P(*)  ■  i  ^  log  (z  -  z^)  -  i  r  (log  r  -f  i  ^  ) 

2Tr  2  ff 

for  z-z  e  re^^ 
o 

where  ^  •  -  JLeis  the  velocity  potential 

2  rr 

Y  ■  .  P  log  r  the  stream  function 

2Tr 

A  permanent  flow  limited  by  a  free  surface  near  which  is  located  a 
vortex,  will  be  defined  by  its  complex  potential  whose  stream  function 
will  be  written  : 


(2.1.) 

(2.2.) 

(2.3.) 


Y  •  Voy  r  log  r  ^ 

2r» 


Voy  +  y 


(2.4.) 
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The  ^  function  is  called  perturbation  stream  function.  It  should 
fulfill  the  boundary  condition  ex-ressed  in  (l.15.).In  this  purpose, 
ve  shall  consider  the  4^  function  as  a  result  of  the  complex  poten¬ 
tials  generated  on  one  hand  by  the  vortex  itself,  on  the  other  h  ad 
by  the  image  vortex  and  finally  by  a  complementary  function  ^  . 

We  can  thus  write  ; 


=  P  log  r  -  1  log  r^  +  ^ 

2V  2  V 


(2.5.) 


^  =  P  log  r  +  y. 

aTT 


(2.6.) 


with  r  =  +  (y  +  f)^  and  r^  =  \/^  +  (y  -  f)^  (2.7.) 

/oo  . 

We  suppose  that  A  ~  j  ^  ("Oe  ^  cos  Icxdk  (2.8.) 

'  /o 


(2.8.) 


where  04  (k)  is  a  function  of  k  to  determine.  Starting  from  the 
expressions  (2.6.),  (2.8.)  and  from  the  condition  (l,15.),  we  obtain 


f  04.  (k)k  e^^  cos  k  x  d  k  +  P  ( 
'  0 


2ir 


+_£  -  y  -  f) 

2  2 


(2.9.) 


[  j  OC  (k)e^^  cos  kx  d  k  P  '-'g  1 

^  J°  o  TT  „  " 


Now,  at  the  free  surface,  y  =  0  and  =  r.  The  expression  (2.9.) 
becomes  : 


/  00  /oo 

y  Ot  (k)  cosk  X  k  d  k  +  T  _f_  »  __g_  (  (X  (k)  cos  k  x  d  k 

IT  2  „  2  y 

*'  Vo  /o  \ 

(2.10.) 
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and  admitting  that  : 


SL  j.  - 

TT 


e  cos  kx  d  k 


(2.11.) 


we  obtain  : 


I  C(  (k)  cos  k  X  k  +  _^  e  cos  k  X  -  _ja_  •<(k)  coskx  jdk  =  0 

0  IT  Vo^  (2.12.) 


Thence  the  function  (k)  : 


«  (k)  = 


r  e  -  kf 

V  (I  -  k) 


wi  th 


g  -  g 

Vo^ 


(2.13.) 


If  this  expression  of  0^  (k)  is  inserted  in  (2.6.),  the  perturbation 


stream  function  can  be  written  ; 


JLlogr.J:  /® 

>17  r  V 


e  >  kf  e^^  cos  k  X  d  k  (2.14.) 
(r  -  k) 


According  to  the  sliding  condition,  the  free  surface  is  defined  by  : 


^  =  -  JL 


(2.15.) 


From  (2.14.)  and  (2.15.)  and  taking  into  account  the  conditions 
previously  set  forth  for  the  free  surface,  we  obtain  ; 


Vo  nv„  (I- 


kf  cos  k  X  d  k 


(2.16.) 
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The  integral  is  undetermined,  but  if  x  is  positive,  its  principal 
value  is  equal  to  the  real  part  of  the  expression  (9) 


.Tl  -:f  +  iKx 
1  e  +1 


00 


e  -  imf  -  m  dm 
(im  -  K) 


(2i17.) 


and  we  have  : 

a  =  -L 

Vo 


-If  .  «  r 

e  sin  Kx  -  / 


Vo 


IT  ^ 


(g  sin  mf  -  m  cos  mf)e 


-mx 


dm 


2  2 
r  +  m 


(2.18.) 


The  function  (2.16.)  is  an  even  function  in  x.  For  x  ^  0,  we  shall 
have  : 


(g  sinmf  -  m  cos  mf)e*"^  dm 

2  ^2 
m  +  g 


(2.181) 


If  we  add  a  system  of  stationnary  waves  to  this  solution  ; 

l^a  =  -  r  sin  gx  (2.19.) 

Vo 

in  order  that,  in  accordance  with  the  physical  reality,  the  free 
surface  at  upstream  of  the  profile  has  no  wave,  the  expressions 
(2.18.)  and  (2.18')  become  : 


r' 


X  Lo 

which  give  the  solution  to  be  obtained. 


2,3t  Free  surface  for  an  immerged  source . 


The  logical  process  which  leads  to  the  determination  of  the  free 
surface  equations  for  an  immerged  source  is  quite  similar  to  all 
that  has  been  exposed  for  a  vortex. 


The  flow  around  a  source  in  the  presence  of  the  free  surface  is 
defined  by  its  stream  function  : 


*4^  =  Vo  y  +  Q  0  +  X  =Voy+4^ 

2  Ti 


(2.21 .) 


The  perturbation  stream  function  must  fulfill  the  condition  (l.t5.). 
We  suppose  that  ^  is  the  total  of  the  potentials  resulting  on  one 
hand  from  the  source  itself,  on  the  other  hand  from  its  image,  and 
finally  from  a  complementary  function  ^  , 


We  suppose  that 


foo 

X  =  y  (k) 


ky 

e  sin  kx  d  k 


(2.23.) 


where  cX.  (k)  is  a  function  to  determine. 


The  perturbation  stream  function  can  then  be  written  : 

'f  =  Q  (0^  +  ©  2^  Jo  ^  ^  ^  ^  ^  (2.24.) 

2  ^ 

The  first  term  of  the  second  member  is  the  imaginary  part  of  the  complex 
function. 
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Thence  : 


J _ ^  (0^  +  Og) 

ay  2T1 


R  _1 _  +  1 

27)  z  +  f  z  -  f 


R  _2_  cos^?  ^  -isinO  ^  +  cos  Sg  "  isin  0^ 

2ir 


(2.26.) 


At  the  free  surface,  y  =  0  and  0  .  =  b’  ;T.  =  T’-,  we  have 

1  2  1  2’ 

^  Q  (  k  X  d  k  (2.27.) 

a  y  2  I*  /)  r^  ^  ° 


(2.27.) 


Taking  account  of  the  conditions  expressed  hereabove,  we  shall  replace 
(2.24.)  and  (2,27.)  in  (1.15.)  : 


2  J  sin  kxdk  ■¥  J  oi  (k)ksin  kx  :  -  i  o(  (k)sin  kxdk 


(2.28.) 


from  which  we  can  obtain  the  function  oi  (k) 


«  (k)  =  -  5  e  -  kf 

Ti  k  -  K 


where  K  = 


(2.29.) 


At  the  free  surface,  th^  ,xpr»;ssion  (2.24.)  will  bv.  written  : 


/  _i  e  -  kf  sin  kxdk 

IT  (k  -  r) 


(2.30.) 
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The  principal  value  of  the  integral  is  equal  to  the  imaginary  part 
of  (2.17.).  We  shall  thus  obtain  : 


cos  Kx  -  _2_  / 

Vo  Vo  fl  '  ^ 


(leos  mf  -f  m  sin  mf)  e"'"*  am  (2.31.) 


2  ,2 

m  +  K 


for  ^  >  0 


is  an  odd  function  in  x,  for  x  ^  0,  we  shall  write  : 


iM  -Iff 

=  __Q  e  '^"'cosKx  +  Q 


(K  cosmf  +  m  sinmf)  dm  (2.31.') 


2  2 
m  +  r 


To  cancel  the  waves  at  upstream  infinity,  we  must  add  a  series  of 
waves  of  the  following  shape  : 


^  -Kf 

Q  e  cos  K  X 


(2.32.) 


It  comes  for  x  >  0  and  x  Z.0  I 


h  =  ^  e”^^cos  Kx 
^  Vo 


rivo  /o 


(K  cos  mf  +msin  mf)  e””**  dm 
- - 


00  y 

(K  cos  mf  +  m  sin  mf)  e  dm 

2  ^2 
m  +  K 


(2.33.) 


2.3.1.  Starting  from  the  expressions  we  have  just  determined,  we  shall 
obtain  the  expression  of  the  free  surface,  due  to  an  immerged  doublet  : 


n  D  •  1 


source 


at  point  ^x+  ^  sinks  at  point  0  =  d  ^  yi 


^  if 

D  *  -  Q  d  %e  ^^K  sin  Kx  -t-  Q  A  x  j  (mK  cosmf  +  m^sinmf)e  ^^dm 

r.  /o  2  2 


II  'o 


2  _2 
m  +  K 


(2.34.) 


.  ^ 
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- 

f 

4 

1 

8 

20 

-  5,0 

+  0,0700 

+  0,1198 

’  4,5 

+  0,0800 

+  0,1448 

-  4,0 

+  0,0900 
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-  3,5 

+  0,0500 
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-  2,5 

+  0,0830 
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-  2,0 
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-  1,5 

+  0,1250 

+  0,2110 

+  0,3676 

-  1,0 
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+  0,2670 

+  0,4453 

•  0,5 

+  0,1750 

+  0,3560 

+  0,5649 

-  0,4 

+  0,1820 

+  0,3750 

+  0,5978 

-  0,3 

+  0,1880 

+  0,3970 

+  0,6286 

-  0,2 

+  0,1960 

+  0,4180 

+  0,6627 

-  0,1 

+  0,2040 

+  0,4380 

+  0,6965 

'v,00 

+  0,2080 

+  0,4520 

+  0.7305 

0,1 

+  0,2090 

+  0,4630 

+  0,7638 

0,2 

+  0,2000 

+  0,4740 

+  0,7959 

0,3 

+  0,1820 

+  0,4810 

+  0,8259 

0,4 

+  0,1540 

+  0,4850 

+  0,8517 

0,5 

+  0,1190 

+  0,':780 

+  0,8781 

1,0 

-  0,150 

+  0,3720 

+  0,9441 

1,5 

-  0,419 

+  0,1350 

+  0,9341 

2,0 

-  0,518 

-  0,170 
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And,  considering  that  tkese  complex  potentials  can  be  written  in 
the  following  way  : 


0  AX 

=  Vo  b^ 

(2.35.) 

271  z 

z 

2  7r  » 

0  ^  X 

(2.36.) 

Vo 

and  then  : 


/  e  sin  m  f  d  m 


(2.37*) 


The  expression  (2.34*)  becomes  : 
D  a  2  b^  Ke"*^^sin  lx  +  2b 


+  2  b 


^  Isl-UL^  mfe'”^  dm  + 

y°  +  K'' 

(  -  r^sin  mf  tm  co3mf)e~'"^ 


2  2 

m  +  ir 


(2.33.) 


where 


P  =  2b^  f  -  2l?b^  re"^^sin  Kx  -  2b^I  j  (K  sinref  -  K  mcosmf )e~”*^dm 


for  ^  ^  0  and  V  /  0 
/ 00 

►\  =qAX  /  ^  niir  ,'nc  1 


(2.39.) 


Vo  II 
*2  / 


/  00 

'  ^  /  (itiK  cos  mf  ■*■  sin  mf)e"*^  dm  • 

Ti  ml  -I.  KV 

I  ^  2  2  2  ^  2 

Jq  e  sin  mf  +  2b  {-TTsLnmt  +  Km^o^mf) 

i,  -  2b^t  f  (Ksin  mf  -  mcos  mf)e”*^m  (2.39 

'''•  /o  _2  .  ,.2 


tT*l  +  Ifl 


m/, 

e  dm 


2  2 
*  2b‘^  f  -  2b^i: 


(2.39') 


m  +  K 


Both  equations  are  similar  to  those  given  by  Lamb  (9)  for  the  case  of 
an  immerged  cylinder. 
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2,4.  Numerical  Calculations. 


The  main  difficulty  to  apply  the  formulae  we  have  just  established, 
consists  in  calculating  the  integral  terms  of  formulae  (2;20.)  and 
(2.33.)  .  With  an  arithmetic  computer  ,  we  have  determined  the 
values  of  these  integrales  for  various  values  of  the  parameter  X  i 

f 

Prom  these  values,  we  have  easily  computed  the  shape  of  the  free 
surface  of  a  flow  perturbed  by  a  source  or  a  vortex  for  various 
values  of  the  parameter  f  .  The  tables  2.1.  and  2.2.  and  the  diagrauns 

A 

2.1.  and  2.2.  give  the  numerical  values  and  the  aspect  of  these 
surfaces . 


®  The  calculations  have  been  effected  by  the  Laboratoire  de  Calcul 
Num^riquc  de  I'Institut  Blaise  Pascal. 
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3.  Analog  Determination  of  the  camber  effect  of  the  mean  line- of  an 
inunerged  profile  submitted  to  an  imposed  pressure > 

3.1.  Definition. 


The  profile  having  a  chord  S  and  an  infinite  span  is  represented  by 
its  projection  on  the  x's  axis  at  a  depth  y  =  -  f  starting  from  the 
static  level  of  the  free  surface  coinciding  with  the  x's  axis.  The 
y's  axis  divides  the  profile  into  two  equal  parts. 

The  free  surface  of  the  flow  will  be  characterized  by  its  wave  length  A 
at  downstream  infinity.  The  flow  will  be  characterized  by  the  ratios 
A  and  f  .  The  Froude  number,  with  respect  to  the  depth  or  to  the 

f  * 

chord,  can  be  deducted  as  a  function  of  ^  and  f  : 

f  S 

^  f  =  2jl  f  (3.1.) 

Vo^  A 

Vo  =  \/  _1_  A  (3.2.) 

V  gf  aTT  f 

where  K  is  the  argument  of  the  trigonometric  function  in  (2.20.)  and 
(2,33.),  and  Vo  the  velocity  of  the  uniform  flow. 


and 


Kf  = 


Pf 


3.2,  Set-up  of  the  problem  and  chosen  pressure  law  ; 

Dealing  with  the  inverse  problem  permits  the  determination  of  the 

shape  of  the  mean  line  of  a  profile  for  a  lower-upper  surface  pressure 

difference  distributed  on  the  chord,  so  as  to  avoid  the  infinite 

overvelocities  at  the  leading  edge.  We  know  that,  in  this  case,  the 

chord  of  the  mean  line  is  placed  at  a  given  ideal  angle  of  attack. 

The  lift  coefficient  Cj^  of  the  profile  will  therefore  be  the 

"adaptation"  C, , 

L 
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The  adopted  load  coefficient  distribution  is  that  of  the  mean  line 
NACA  65  We  can  write  : 


A  Cp  =  C'  .  t  S  =  C 

*■  r  ‘‘IT 


■  ^  \/l  -  4  (  X)' 
t  IT  V  s 


2  ^  2 
(3.3.) 

where  ^  represents  the  vortex  density  along  the  chord..  In  a  A 
interval,  the  vortex  density  Jf"  makes  appear  a  vortex  which  we  shall 
call  AT  ,  the  value  of  which  is  : 

K 


s  s 


In  an  indefinite  fluid,  the  mean  line  corresponding  to  this  vortex 
distribution  is  given  by  the  expression  ; 

1*  L.  "  -1.  1  -  4(i)^  -iLtLi.  (3.5.) 

S  O  2  IT  S  2  2 

w 

Figure  3.1.  gives  us  the  shape  of  functions  (3.4.)  aoid  (3.5.). 


3.3.  Boundary  conditions  and  formulae. 


a)  On  the  segment  representing  the  profile,  the  boundary  condition 
belongs  to  the  Neumann  type .  Replacing  A  by 

3  y  dy 

and  talcing  into  account  the  sign  of  the.  normal,  it  can  be  written  : 

-  .r  i  /  -  ./C^\2 


n 


iJi  \/  1  - 

TTs  s 


(3.6.) 
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b)  When  replacing  the  profile  by  N  vortices  located  along  its  projec¬ 
tion,  the  geometrical  shape  of  the  free  surface  will  be  considered  as 
the  superposition  of  the  effect  of  each  punctual  vortex  ;  the  contri¬ 
bution  due  to  the  vortex  of  ^  P intensity  and  located  at  point 
X  =  t  n,  y  =  -  f  is  K  ,  jf  n, 

Vo 


According  to(2.20),  we  have  : 

In  =  -  2e"^'^sinK(x  -T  )-1  /  (isinmf  -  mcosmf)e”"'^  ^ 

"  fi/o  ^ 

>«  >  0 


)f  =  -1 

.?n) 


The  free  surface  will  therefore  be  ; 
n  =  N 

h  ,  V  “  2  *1  (x,  (F,  n)  ^  ^  t  n  (3.8.) 

7  (?»  )  n  =  1  "  „ 


And,  talcing  into  account  (l.12.) 

n  =  N 

=  -  Vo  ^  -  I  ^  (x,  AJlia  (3.9.) 

P  p  n  »  1  p 


representing  the  values  of  the  potentials  to  be  set  at  the 
free  surface . 


c)  For  X  «  -  oo ,  the  perturbation  potential  is  equal  to  zero,  therefore  : 


I'OO 


(k  sinmf  -  m  cos  nif)e 


">(«  -?„) 


(3.7.) 


m 


<0 


V  (-CX3  ,  y)  =  0 
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d)  For  a  dot  rather  far  from  the  profile,  the  only  term  which  remains 
in  the  expression  (3.7.)  is  : 

~  ^  n^* 

For  the  whole  symmetrical  lot  of  vortices,  we  obtain,  after  operating 

n  =  N 

=  2  e”*^^  sin  Kx  ^  ^  ^  ^  ^  n  .  cos  K  ^ 

r  n  =  1  r 

If  we  take  a  maximum  or  a  minimum  of  this  wave,  i.e.  an  at"  .issa 

X  =  (2n  +  l)  ^  ,  the  vertical  component  of  the  perturbation  velocity 
4 

is  equal  to  zero,  therefore  : 


_ Lf,..a-x}  =  0 

e)  For  y  =  -  00 ,  we  suppose  that  the  flow  is  not  perturbed  ;  we 
shall  have  : 


=  0 


3.3.1 

S'! 


.  We 
and 


shall  consider  the  function  as  the  sum  of  two  functions 
'f  2  determined  by  the  following  boundary  conditions  : 


'‘-•If./? 
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Func'jior.  S'"  , 


a)  At  'the  free  surface 

.,(x,0)  -  0 

b)  Cn  tbe  profile 

Ti  o  b 

c)  for  X  =  -  CD 

S’^  .,(-00  ,  y)  =  0 

d)  for  X  =  CO 

''-y  -jioo,  y)  =  0 


e)  for  y  ■■■  -  00 
I  .,(x,Co)  a  0 


a') 

WjC*,  0)  .  I 


n  =  N 


n  =  0 


b') 


y  ^  c 


2(-  »  y)  =  0 


d')  For  X  "  ( 2n  +  1 ) 'i  :i  -  0,1, 

■f, 

.0 


e' ) 


vy  u.”' 

Th.e  i  " ;  i-urcmcr-'L  of  the  Dotent'  lo  i  ^  and  1  „  on  t’le  electrodcr 

1  2 

cf  'chc  yl- ca  allovs  to  deterir.ine,  as  we  shall  sec  further  't,  the 

shape  of  the  neari  line  of  tlie  profile  'or  each  pair  of  parainatcrs 

and  _f  . 

A 

Fig.  3.2.  sho-.'.'s  the  boundary  liirdts  exposed  hereabove. 
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3.3.2,  Physical  interpretation  of  .y,  isa-ys- 

The  ^  ^  function  superposed  to  a  uniform  flow  represents  a  flow 
perturbed  by  an  immerged  profile  (either  a  mean  line  or  a  non  lifting  thicl 
profile),  the  free  surface  of  which,  according  to  the  assumption, 
remains  rigid.  As  a  consequence,  the  perturbation  function  ^  ^  at 
the  free  surface  is  equal  to  zero. 

On  the  other  hand,  the  g  function  superposed  to  a  uniform  flow 
represents  a  flow  from  which  the  profile  has  been  taken  off,  but  whose 
free  surface  shows  the  same  shape  as  if  the  profile  was  still  in 
the  flow. 

The  image  of  the  total  flow  is  obtained  by  superimposing  the  fields 
of  the  stream  lines  created  by  the  two  above  mentioned  flows. 

Figure  3.3.  represents  the  three  quoted  flows. 


3.4,  Analog  representation  and  realization  of  the  experiments. 

In  order  to  represent  an  important  part  of  the  study  field  of  ^  , 

the  electrical  potential  representing  'f  is  studied  in  a  rheoelectri- 
cal  set-up  made  up  of  an  ordinary  tauik  (tank  z)  and  of  am  other  tamk 
(tank  Z)  representing, in  a  mathematical  inversion,  the  outside  part 
of  the  first  one  (11).  The  two  tanks  are  limited  by  cylindrical 
surfaces  paved  with  electrodes,  connected  by  pairs  one  to  the  other, 
so  as  to  establish  the  same  potential  on  each  homologous  dot  of 
the  two  surfaces.  The  power  of  the  inversion  transformation  is 
given  by  the  product  of  the  radius  of  each  circle  : 


2 

P  ■  oa  X  o'a'  ■  292$  cm  • 
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At  the  axis  center  of  the  tank  Z,  we  put  half  a  clover  of  electrodes 
made  up  of  three  cylinders  having  a  45  mm  diameter  and  tange:  t  one 
to  the  ?,  Each  branch  represents  one  of  the  physical  limits 

of  the  field.  This  set-up  permits  the  representation  of  a  region 
included  between  -  650  cm  x  /^  +  650  cmi 

The  profile  is  made  up  of  a  plate  of  plexiglass  (3  mm  thick,  and 
20  cm  long)  bearing  1  cm  wide  electrodes  on  both  faces.  The  electrodes 
are  conrected  two  by  two  so  as  to  fulfill  the  condition  : 


r 


+' 

n 


Ve  shall  therefore  obtain  : 


(3.10.) 


(3.11 .) 


The  feeding  of  the  electrodes  of  the  plate  is  effected  by  the  large 
resistance  method  (12)  utilized  to  represent  the  Neumann  condition. 


The  free  surface  of  the  flow  is  represented  in  the  tank  z  by  a 
plexi  ■’.?.ss  plate  paved  all  along  with  2  cm  electrodes.  In  the 
tank  Z,  the  part  of  the  surface  representing  the  flow  at  the 
upstream  part  of  the  hydrofoil  is  covered  with  a  continuous  electrod. 
When  necessary,  it  is  possible  to  take  this  electrod  off,  and  to 
carry  on  the  experiment  with  discontinuous  electrodes.  The  downstreaim 
part  is  paved  with  1,5  cm  electrodes.  The  potential  setting  at  the 
free  surface  is  effected  by  the  mean  of  a  couple  of  ne  thousand 
Xf.'es  divider  bridges  mounted  symmetrically.  The  zero  value  is 
obtained  by  the  mean  of  a  transformer  T  permitting  to  inject, 
between  the  t’-o  zero  marks  of  the  one  thousand  tapes  divider 
bridges,  an  electrical  current  having  a  direction  opposite  to 
the  current  direction  through  thj  leads  connecting  the  bridges. 

This  intensity  is  regulated  by  a  potentiometer  mounted  as  a  rheostat. 
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3.4.1.  Representation  of  the  function  ^  ("'Lgure  3.4.) 

In  the  displacement,  the  flow  is  not  p-,  rturbed  upstream  and  downstream, 
condition  3.3.1.  c  and  d.  Therefore,  the  arcs  of  circles  d'c*  and  f'c', 
figure  3.4.,  must  be  covered  with  a  continuous  electrod  connected  to 
the  zero  of  the  measurement  bridge.  At  the  free  surface,  the  imposed 
condition  is  3.3.1.  a,  the  electrodes  aOb,  a'd'  and  c'b'  are  also 
given  a  zero  value.  At  a  depth  important  enough,  the  condition  3.3.1. 
is  fulfilled  by  a  continuous  electrod  on  e'f  connected  to  zero. 


The  condition  3.3.1.  b  on  the  profile  is  fulfilled  by  the  mean  of  a 
large  resistance  set-up.  Let  us  consider  (fig.  3. 4. a)  two  electrodes 
connected  one  to  the  other  and  fed  by  a  potential  -  100  through  a 
resistance  having  a  R  value.  The  intensity  given  by  the  electrodes 
will  be  : 

1  I  =  1  Vin  -  ^ 

<r  ^  h  «r£  h  R 

«  ^  f(  txi)  (3.12.) 

S  S 


=  conductivity  of  the  water  ^  cm”^) 

=  width  of  the  electrodes  (cm) 

=  height  of  water  (cm). 

If  y  is  determined  as  above  mentioned,  to  the  value  -  100  and  if 

we  give  the  resistances  R  rather  high  values,  it  is  always  possible 

to  measure  a  potential  included  between  0  and  -  5.  Now,  let  us  make 

R  inversely  proportional  to  the  values  of  f( _ ),  then  ; 

S 


where  fT 

i 

h 


1  R  •  R 

f{]Ln) 

S 


(3.13.) 
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Tbi  condition  (3«12.)vill  be  written,  if  ve  suppose 

V  -  ~  100  ,  (3il4.) 

and,  if  we  chaoigc  R  by  its  value  obtained  in  (3.13.)  : 

w 

1  f(  S  )  100  »  £.  f(  C  n)  (3.15.) 

0^  h  K  S  S 

This  relation,  where  all  the  participating  factors  are  known, 
permits  to  compute  the  value  of 


r  «  100S  (3.16.) 

^  hj 

The  value  of  K  is  obtained  by  making  several  step  by  step  tests  untill 
we  verify  that  0/_  4^  -  5.  Practically,  two  tests  are  enough. 

If  we  measure  the  potentials  on  each  electrod  of  the  plate,  ve  can 
write,  according  to(l.12.) 


W  -  =  -  Voy  ,  (3.17.) 

^1^'n  '1A  Qn 

p 

To  obtain  a  non-dimensional  result,  we  divide  this  expression  by* 

f  1  n  -  S'  1A  °  >'lVn  (3.18.) 

r  r  c, 

and  in  multiplying  numerator  and  denominator  of  the  second  memb'  r 
by  the  chord  S,  and  then  in  talking  into  account  the  relation  : 


VoS  »  2 


(3.19.) 
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Ve  get  : 


1  n  =  -VlA  =  y1  n  .  2  (3.20.) 

r  r  ® 

3.4.2.  Representation  of  the  function  ^ 

In  the  displacement  ^  the  electrodes  at  the  free  surface  on 
aob,  a'd'  and  c'b'  are  connected  to  the  potentials  calculated  accor¬ 
ding  to  (3.9.)  (condition  3.3.1.  a').  At  upstream  infinity,  the 
flow  is  not  perturbed  (condition  3.3.1.  c')  we  can  thus  impose  the 
zero  potential  on  d'e'  as  well  as  on  e'd',  under  the  condition  that 
the  profile  be  far  from  this  region.  At  a  rather  important  depth, 
the  condition  3.3.1.  is  valid  and  we  impose  the  zero  potential  on 
e'f.  Finally,  on  f'e',  according  to  the  condition  3.3.1.  d' ,  the 
surface  will  have  to  be  insulating. 

On  the  plate,  the  condition  3.3.1.  b'  is  fulfilled  if  the  electrodes 
are  not  fed.  On  these  electrodes,  we  shall  measure  an  electrical  po- 
tential  proportional  to  ■  2  C  n  and  we  can  write  : 

at  n.  Vstn-  ^  ZA.-yg  :£  n  .  2_ 

r  p  ® 

(3.21.) 


2b 


3.4.3.  Composition  of  ^  ^  and  ^ 

The  dimensionless  variables  obtained  by  the  formulae  (3.20.)  and 
(3.21.)  give,  to  an  additional  constant,  the  value  of  the  mean  line 
data  at  the  abscissa  ^  ^ ,  for  each  of  the  displacements  ^  ^  and 
^  g.  The  total  of  the  expressions  (3.20.)  and  (3.21.)  allows  to 
obtain  the  shape  of  the  mean  line  in  a  flow  equal  to  the  sum  of 
and  ^  g.  Therefore  : 

p  r 


V  ' 

11. 

6  n  +  ?'2  If  n 

=  -  2  y  ^  n  (3.22.) 

. 

S  S  . 

3.5.  Experiments  and  results . 


We  have  realized  two  series  of  experiments  for  two  different  ratios 
of  f/s .  These  ratios  have  been  chosen  for  the  following  reasons  ; 
for  depths  superior  to  a  chord,  the  influence  of  the  free  surface 
seems  rather  low,  while  if  the  depth  decreases  by  more  than  half  a 
chord,  we  are  more  and  more  leaving  the  practical  cases.  For  each 
ratio  of  f/s,  we  have  realized  seven  experiments  for  different  £//\ 
ratios  ;  oo  ,  l/4,  l/8,  l/20,  l/50,  l/lOO  and  0, 

The  two  extreme  cases  oo  and  0  are  of  interest  as  they  correspond 
to  two  flows  whose  physical  characteristics  are  also  extreme  ; 
infinitely  heavy  and  non -heavy  fluid. 

Fig,  3.5.  shows  the  shape  of  the  mean  line  which  bears  a  pressure 
law  imposed  for  those  two  boundary  cases  and  for  various  ratios  of 
f/s.  It  is  of  interest  to  point  out  the  lift  coefficient  variation 
of  a  NACA  65  mean  line,  the  angle  of  incidence  of  which  is  equal  to 
zero,  with  respect  to  the  lift  coefficient  in  an  indefinite  fluid. 
Fig.  3.6.  shows  the  unfavorable  effect  of  a  non  rigid  free  surface 
on  the  lift.  Fig.  3.7.  and  3.8.  show  the  aspect  of  the  mean  lines 
obtained  through  the  analog  computation,  for  various  ratios  of  f  . 
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4*  Analog  Determination  of  the  Shape  of  a  Thick  non  lifting  Profile 


when  the  thickness  lav  is  imposed. 


4.1 .  Bet-up  of  the  problem  and  chosen  pressure  law. 


Which  has  been  defined  as  thickness  effect  in  the  illimited  two- 
dimensional  flow,  now  loses  its  simplicity  and  its  meaning  since 
the  symmetrical  thick  profile  constituting  the  hydrofoil  near  the 
free  surface  is  exposed  to  a  velocity  dissymmetry  between  its 
upper  surface  and  its  lower  surface,  even  when  its  angle  of  inci¬ 
dence  is  equal  to  zero,  Obviously,  this  dissymetry  depends  on 
the  shape  of  the  free  surface  defined  for  each  pair  of  the  values 
f  and  f  . 

/\  s 


We  intend  to  find  out  the  deformation  to  be  given  to  a  symmetrical 
thick  profile  in  order  to  avoid  th-.  lifting  effect  °  which  appears 
when  it  is  located  in  a  flow  limited  by  the  free  surface. 


We  have  chosen  as  thickness  law,  because  of  the  assumption  of  the 
small  perturbations,  the  thickness  law  which  gives  a  sharp  profile 
at  both  ends.  The  infinite  slopes  at  the  leading  edge  are  thus  avoided. 
Besides,  the  thickness  distribution  is  symmetrical  with  respect  to 


the  y's  axis.  It  can  be  written  as  follows  : 


=  1  -  4(^)‘ 

S 


- i  i 

o  ^  o 


(4.1.) 


®  In  fact,  the  symmetrical  profile  at  zero  incidence  located  near 
the  free  surface  of  the  flow  will  have  such  a  pressure  distribution 
that  the  lift  will  be  directed  downside. 
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where  s  represents  the  chord  of  the  profile  ;  the  origin  of  the 
abscissae  is  taken  at  the  middle  of  the  profile* 

4.2.  Boundary  conditions  and  formulae t 

a)  On  the  segment  representing  the  profile,  the  perturbation  strecim 
function  takes  on  values  determined  by  the  chosen  thickness  law» 
Therefore  : 

.a'  _\s)  ~  .  + 

^  =  y  -  y  =  e  (4.2.) 

Vo  Vo 

and,  if  we  consider  the  ratios  of  these  values  to  the  maximum 
thickness  : 

1  Ay  -  _ 1_(  7  *  -  'f  -)  -  1  -  4(1)  -S^'g^S 

'  S  %  2 

Vo  e  max  Vo  e  max  ^ 

(4.3. ) 


b)  By  replacing  the  profile  by  N  sources  and  N  sinks  set  along 
its  projection  and  symmetrically  with  respect  to  y,  the  wave  shape 
of  the  free  surface  will  be  considered  as  the  superposition  of  the 
effect  of  each  punctual  source  and  sink.  The  contribution  due  to 
a  source  or  a  sink  having  a  mass  flow  equal  to  A  Q  and  located 
at  the  dot  ;  x  =  ^  n,  y  =  -  f ,  is  1*^  ^  ^  . 

According  to  (2.33.),  we  have  : 


^00 

jo  /  y  =  h  Vo  =  2e  cosK(x-^n)-1_  /  (Kcosmf  +  msinmf)e 

C  T\  - 2 - 2 - ^ 

y>  0  " 

✓ 

(Kcos  mf  sin  mf)e'^^  ^  ”  ^”^dm 


-m(x- 


to  =f7  12  =  1 
‘/JQ  ff 

X  Z_o 


m2  +  K 


(4.4.) 
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where  /S  Q  represents  the  mass  flow  of  a  source  equivalent  to  the 
profile  (see  §  1.3.)  between  any  two  ordinates.  We  shall  define 
the  value  of  this  mass  flow  as  a  function  of  (4.1.). 


e 


e 


max  4 

S" 


+  — 

y  -  y 

4.5. 


The  values  of  the  perturbation  velocities  v"^  and  v  can  be  obtained 
when  taking  into  account  the  sliding  condition  (l.8.),  and  in 
derivating  (4.5.) 


V  -  V 


^  y~)  =  -  ^  Vo  e  max  ^  (4.6) 

3^ 


The  mass  flow  of  an  elementary  source  will  be  given  by  ; 


dQ  =  (v^-v’)d^  (4.7.) 

and,  by  integrating  between  ^  ^  ^  and  ^  ^ 


+  iAt 


A 


-  8  Voe  max|P  d  1^  =  -8  Voe  max  ^  ^  ^  ^ 


(4,8.) 


5 

If  we  replace  (4.8.)  in  (4.4.)  to  estimate  the  wave  motion  of  the 
free  surface,  and  if  we  take  into  account  (l.12.),  we  obtain  a 
Dirichlet  condition  for  the  free  surface  : 

(X)  -7^.  dn)  Aalj!  "  -IlYv 


Vo 


(4.9.) 
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We  can  write,  to  make  this  expression  dimensionless 


It  tY!  0>f 


(4.10.) 


and  if  we  make  the  sum  for  all  the  sources  and  sinks  which  replace 
the  profile  ; 


n  =  2N 


1  rum  .  -1,  A? 


n  =  1 


e  max 


n=N 


n=2N 


^  *7  (x,  f  2-  ^''(^(x,^n)]fnj 

n=1  '  n=N+1 


(4.11.) 


If  we  take  into  account  the  fact  that  for  the  sources  and  sinks  the 
absolute  value  of  )r  is  the  same  but  of  opposite  sign,  we  dial 1  have, 

w  n 

for  a  symmetrical  distribution  of  sources  and  signs  : 


n=N 


^  -  !2iiL)= 

,2 


Voe  max  n=1  e  max  S 


n=N 


n=2N 


X»V^(x,!rn)  "21  n’^(x,  -J?  ) 

n=1  n=N+1  ‘‘ 


(4.12.) 


The  first  member  represents,  to  a  constant,  the  values  of  the  potential 
to  be  set  at  the  free  surface. 


c)  for  X  =  -  00 ,  the  perturbation  potential  is  equal  to  zero,  therefore  : 

(-  oo  ,  y)  =  0 

d)  Like  in  3.3.d)  at  downstream  infinity,  the  only  term  which  remains 
in  the  first  expression  of  (4.4.)  is  : 


=  "^^(K  ,']fn)  ^  ^  k(x  -  t.  n) 


Vo 


(4.13.) 
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For  the  whole  symmetrical  lot  of  sources  and  sinks,  we  obtain  ; 

^  V  2e~^^  cosK(x-*Jn)-cosi:(x+  J^n) 

*1  (  >  ,  JT  n)  ^  ^ 


=  4e~^^sinK  n  sin  Kx 

And  finally,  according  to (4. 12.)  : 

n  =  2N 

t  .  32  e  -“a?  V  t  n  sinKx  sinK '|£  n  (4.14. 

Voe  max  n  =  N  +  1 

The  maximum  or  minimum  value  will  be  obtained  for  x  =  (2n+l)  A  . 

4 

From  this  abscissa,  the  stream  function  will  be  so  that  ; 

liL  =  0 

J  X 

as  the  vertical  component  of  the  perturbation  velocity  is  equal  to 
zero . 

e)  for  y  =  -  00  we  suppose  that  the  flow  will  not  be  perturbed  ; 
the^function  will  therefore  be  equal  to  zero. 


«  0 
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4.2.1  .  We 
‘K  1  and 


shall  consider  the  'y  function  as  the  sum  of  two  functions 
^  5  which  are  ruled  by  the  following  boundary  conditions  ; 


Function  ^  ^ 
a)  at  the  free  surface 


b)  on  the  profile 

^  1  ~  V  1  =  1  - 

Voe  max  S 


Function  ^  ^ 
a') 

2  (formula  4.13.) 
Voemax 

b') 

;  ..  0 

Voe  max 


c)  at  downstream  infinity 


I 


1 


«  0 


d)  at  upstream  infinity 


c '  )  at  .V  =s  ( 2n+1 

4 

n  »  0,1 ,2. . .n 


d') 


=  0 


X 


1 


=  0 


0 


e)  for  y  =  -  00 


e*) 


1 


=  0 


0 


The  ^  ^  and  ^  2  potential  measurement  on  the  electrodes  of  the  plate 
allows,  through  the  sliding  condition  on  the  profile,  defined  in 
(§1.14.)  the  determination  of  the  ordinates  y"*"  and  y".  Paragraph  4.3.3 
shows  the  operations  leading  to  this  determination. 

Figure  4.1.  represents  the  boundary  conditions  expressed  above. 


33 


4.3.  Analog  representation  and  experiments . 

The  experiments  have  been  realized  in  a  quasi-indefinite  rheoelectrical 
set-up  whose  characteristics  have  been  given  in  3.4. 

4.3.1  .  Representation  of  the  function  ^  . 

The  condition  b)  of  4.2.1.  on  this  profile  is  obtained  by  feeding  each 
pair  of  electrodes  of  the  plate  through  a  transformer  with  a  separate 
feeding  T  and  having  a  l/3  ratio  ;  its  secondary  contains  a  potentiome¬ 
ter  p.  Obviously,  the  electrodes  are  not  connected  one  to  another. 

The  measurement  of  the  potentials  is  effected  by  the  mean  of  a  measure¬ 
ment  bridge  PM,  fed  by  a  transformer  having  a  l/3  ratio,  so  as  to avoid 
a  current  supply  of  the  electrodes  during  the  conduct  of  the  measure¬ 
ment  (iPig.  4. 2. a). 

The  limits  of  the  field  are  covered  with  a  continuous  electrod  connec¬ 
ted  to  the  zero  of  the  measurement  bridge  in  the  same  way  as  in  3.4.1. 

The  voltage  between  two  electrodes  of  the  plate  will  be  proportional  to 

the  expression  (4.1.).  On  these  electrodes,  we  measure  the  values  of 

electrical  potentials  proportional  to  4^  t  w  and  H''  1  ‘g. 

Vo  e  max  Vo  e  max 

+• 

According  to  (l.14.)  and  calling the  reference  potential  measured 
at  the  leading  edge,  we  can  write  : 

v1  t  =  -  ^  1(g  -Y  1A  and  y1  't  =  -  ^  ^  1A 

e  max  Vo  e  max  e  max  Vo  e  max 

(435.) 

4.3.2.  Representation  of  the  function  ^ 

On  the  electrodes  representing  the  free  surface,  wc  are  setting  the 
electrical  potential  values  proportional  to  (4.13.)  -  condition 
4.2,1.  a'  -  On  the  other  limits  of  the  field,  the  conditions  are 
those  of  3.4.2.  (figure  4.2.b), 
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According  to  cofidition  4.2,1.  b*  : 


+ 

2 


t 


2 


Vo  e  max 


Vo  e  max 


where  Y  ?  Y  2 

The  electrodes  of  the  plate  must  now  be  connected  two  by  two. 
We  shall  measure  on  these  electrodes  an  electrical  potential 
which  will  be  proportional  to  V  2 

Vo  e  max. 


We  can  write  that,  to  a  constant  ; 

r  2^  =  y2  k*  =  -  Y  2 C-  V  2A  =  .  V  2!"  Y  2A 

e  max  emax  Voe  max  Vo  e  max 


where 


f 


+ 

2A  represents  the  electrical  potential  -.t  the  leading  edge. 


4.3.3.  Composition  of  the  functions  ^ 


and  VSy 
I 


2 


The  addition  of  formula-  (4.15.)  and  (4. 16.)  permits  to  define 
the  ordinates  of  the  jrofile  in  the  following  way  : 


y  ifiA. 

e  max  e  max 

and  : 


j _ 

Voe  max 


(4.17.) 


1 


(4. 18.) 


e  max 


e  max 


Vo  e  max 


which  respectively  give  the  upper  amd  lower  surfaces  of  a  non 
lifting  profile  for  each  pair  of  parameters  having  a  flow  f 

S 

and  f  . 

4.4.  Experiments  and  results. 

Figures  4.3,  and  4.4.  show  the  results  obtained  for  two  series  of 
experiments  where  the  shape  of  the  profile  is  given  as  a  function 
of  y  and  .  One  can  notice  that  the  dissymmetry  of  the 
e  max  S 

velocities  due  to  the  presence  of  the  free  surface  makes  the  incidence 
angle  increase  :  this  angle  is  generally  small  with  respect  to  those 
due  to  the  lift  effect  for  usual  values  of  e  max/s  and  of  Cj^, 

We  have  noticed  in  the  course  of  the  experiments,  that  the  profile 
remains  nearly  symmetrical  for  f  limits,  0  and  oo ,  and  for  depths 

A 

superior  to  half  a  chord.  This  fact  can  be  explained  in  the  light 
of  the  small  perturbation  assumptions  which  we  had  set  forth  from 
the  beginning. 
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